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ABSTRACT 

Aims. The aim of this paper is twofold: (a) to determine the reliability of infrared (IR) emission to trace star formation in individual 
star-forming sites of M33, and (b) to outline a new method for testing the distribution function of massive stars in newly formed 
clusters. 

Methods. We select 24 /im IR sources from the Spitzer survey of M33 with Ho- counterparts and show that the IR luminosities have 
a weak dependence on galactocentric radius. The IR and Ho- luminosities are not correlated. Complementing the infrared photometry 
with GALEX-UV data, we estimate the bolometric luminosities to investigate how they are related to the Ho- luminosities. We 
simulate a theoretical diagram for the expected bolometric-to-HQ- luminosity ratio, Lj,„//L// b , of young clusters as a function of the 
cluster luminosity. We then compare the observed Lj, //L//„ ratios with the theoretical predictions. 

Results. In the log(L/, /)-log(Li, //L//„) plane, stellar clusters should be born along a curve that we call the cluster birthline. The 
birthline depends on the stellar initial mass function (IMF) at the high-mass end, but not on the cluster mass function. For an upper 
stellar mass limit of 120 M Q , the birthline is flat for L 6o/ > 3 x 10 39 erg s -1 because all clusters fully sample the IMF. It increases 
toward lower luminosities as the upper end of the IMF becomes incompletely sampled. Aging moves clusters above the birthline. 
The observations of M33 show that young isolated clusters lie close to the theoretical birthline for a wide range of L 6o; . The observed 
L&oi/L/fa ratio increases toward low L&,/ like the theoretical curve, indicating that luminosity is not proportional to Ho- emission for 
low mass clusters. The best fit to the birthline is for a randomly sampled IMF, in which the mass of most massive star in a cluster is 
not strictly limited by the cluster's mass, but can have any value up to the maximum stellar mass with a probability determined by 
the IMF. We also find that the IR luminosity of young stellar clusters in M33 is not proportional to their bolometric luminosity. This 
irregularity could be the result of low and patchy dust abundance. In M33 dust absorbs and re-radiates in the IR only part of the UV 
light from young clusters. 

Key words. Galaxies: individual (M 33) - Galaxies: ISM - Galaxies: star clusters - Stars: mass function 
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Corb elli & Schneidedl997tlHoopes & Walte rbos 20001 ICorbelhl 

20031: lEngargiola et alj|2003b iHever et alj|2004l) . Together with 

' „ , , , c ■ . „ j... . r. c recent works on the metallicity content of the young and old 

« Our knowledge or interstellar conditions that favor the birth or „ . . , ,. , , , 

>, , . 4.1 u j i\/r ii nr u stellar population across the M33 disk, they allow us to trace the 

stars and clusters is mostly based on the Milky Way, but these , , , \- ■ , , • ■ ■ , , , ,• , 

,. ... . , „ . . . ... global formation and evolution history of the nearest blue disk 

studies are ill-suited for examining massive stars, which are an ° . rrr — . . . h ' 



, , , „. .,, , ,. rm. ' , t ,,.„ galaxy (Magrini et al. 2007a, and references therein). 

«: ■ embedded population with a short lifetime. The number of Milky ° ' 

03 Way clusters that one can sample is also limited by extinction | Verley et al.| (|2007|. |2008[) have shown that the recent Spitzer 

and distance uncertainties. These problems are somewhat di- ima 8 es of M33 allow us to stud y the average radial variation 

minished in studies of other galaxies, where imaging with high of star formation properties across the disk, and to locate stel- 

sensitivity and resolution is possible thanks to space observato- lar clusters and individual stars with infrared luminosities as low 

ries such as Spitzer, HST and GALEX. Soon we shall be able as 2 x 10 37 erg s 1 (corresponding to a B2 spectral type star 

to know if the IMF and extinction in star-forming sites varies on the main sequence). M33 therefore gives us a unique oppor- 

along the Hubble sequence, and if ISM perturbations from these tunit y t0 study the properties of many individual star-forming 

sites vary with cluster position and mass. In the mean time, sites, containing a range of populations from single O and B- 

the best galaxies for stud y are in the Local G roup. M33, at type stars t0 y oun 8 stellar clusters and associations. The observ- 

a distance D = 840 kpc (iFreedman etallfll9l . has a higher able l uminosity range extends up to (4.5 ± 1 .5) X 10 7 L H for NGC 

SFR per unit ar ea than M31 (3.4 versus 0.74 M Q Gyr 1 pc 2 , 604 (Israel et al.|U98l, 1 1 986|; [Hippelein et alj|2003|), the second 

iKennicutlJI 19981) and a lower extinction towa rds star-forming re- brightest HII complex in the Local Group. 

gions because of its lower inclination (e.g. Israel & Kennicuti The aim of this paper is to examine the conditions of the in- 

I l980t iMagrini et al.1l2007bl and refer ences therein). Bearing no terstellar medium where stellar clusters are born, to define stel- 

prominent bulge and no signs of recent mergers, M33 is a pro- lar cluster properties related to the birth and death of massive 

totype for star formation in a quiescent galaxy and for evolu- stars, and to examine star formation rate (SFR) diagnostics of 

tionary scenarios in blue, low luminosity systems. Observations individual complexes in M33. We will show the importance of 

of the molecular, atomic and ion ized gas in M33 have been complementing infrared photometry with photometry at optical 

carried out in the past years (e.g. iDeul & van der Hulstlll987t and ultraviolet wavelengths. Previous papers outlined the need 
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for infrar ed photometry to supplement Ha in deter mining the 
SFR (e.g. ICalzetti et allllool llCennicutt et al.ll2007l forM51). 
In bright star-forming regions of more distant and brighter galax- 
ies, such as M5 1 , there is a tight linear relation between IR lumi- 
nosity and the st rength of extinction -corrected, optical recombi- 
nation lines (e.g. lCalzetti et alj|2005l) . On the other hand, in late- 
type galaxies there are strong variations in the Ha-to-infrared 
flux ratio, which limits our ab ility to trace star fo rmation via in- 
frared photometry alone (e.g. Cann on et al. 2006, for the metal 
deficient Local dwarf NGC6822, where the Ha-to-IR flux ra- 
tio variations are as high as a factor 10). UV observations are 
also needed to trace the properties of star formation, particularly 
where the HII regions are density bounded or the O-type stars 
have evolved off the main sequence. GALEX observations of 
M33 have shown that the UV luminosity of a region is compa- 
rable to the IR luminosity (IVerlev et alj|2008l) . and is therefore a 
significant contribution to the total. 

Here we employ a new technique to extract and measure 
sources from M33 images. Because dust and its associated IR 
emission are distributed in both star-forming complexes and 
the diffuse interstellar medium, we do photometry with varying 
apertures to match the source size at each wavelength. We also 
combine IR, UV and Ho- fluxes to give source properties down 
to very faint emission levels. 

It is currently unknown if IR fluxes are good diagnostics for 
the SFR in areas where the bolometric luminosity of the new- 
born cluster is below 10 4() erg s _1 . M33 is an excellent candidate 
for studying individual low-mass clusters because its proximity 
minimizes contamination from the surrounding ISM. Whether 
low-mass clusters contribute substantially to the SFR of a galaxy 
depends on the initial cluster mass function (ICMF). In M33 the 
fragmentation process seems to favor small masses: the mass 
spectrum slope for giant molecular clouds (GMCs) is steeper 
than -2, as is the HII region luminosity function dete r mined 
from Ho- emission (e.g. IWvder et al.lll997l) . IVerlev et ail d2007l) 
have shown that the luminosity function at 24 pm can be de- 
scribed by a double power-law that is shallower than -2 at its 
faint end. While a change of slope like this could mark the tran- 
sition between poor and rich clusters, occurring at the luminosity 
of the brightest star formed, a shallow slope at the low luminos- 
ity end can also indicate a decreasing dust abundance for lower 
masses. Giant complexes may form primarily in dust-rich envi- 
ronments, e.g. close to spiral arms, while low-mass clusters may 
form everywhere, also in low-dust environments. If some low- 
mass clusters are born in a low or patchy dust environment, then 
we cannot use IR luminosities alone to infer the SFR. 

In this paper we will study the large scatter in the IR-to-Ha 
ratio that is observed when young, low-luminosity stellar clus- 
ters are selected in M33. We show that complementing the IR 
photometry with UV photometry helps recover the total cluster 
luminosity and improves the correlation between stellar contin- 
uum and gas recombination lines from star-forming regions. We 
will examine in detail the theoretical and observed relation be- 
tween Lboi and L// a , the deviations from the linear regime, and 
some properties of the massive stellar population in young clus- 
ters. 

We shall also investigate the 8pm emission from PAHs in in- 
dividual HII regions. Generally, PAH emission is low in dwarf 
and low-luminosity blue galaxies compared to other IR radia- 
tion. This may be related to the low metallicity and intense stellar 
radiatio n field in these systems relative to those in spiral galax- 
ies (e.g. iGalliano et alJl2003ULu et alJl2003l:lHouck et al.l b004; 
iHogg et alJ 120051: iRosenberg et all 120061) . iHunter et alj (1200 ll) 
for example, used Infrared Space Observatory (ISO) mid-IR 



imaging and far-IR (FIR) spectroscopy to examine the proper- 
ties of five dwarf irregular systems. They found that PAH emis- 
sion, which is associated only with the brightest H II regions, is 
depressed relative to that of small grains and far-IR. In addition, 
the integrated [C II] emission relative to PAH emission is high 
in dwarfs compared to spiral galaxies, sugg esting that atomic 
carbon is elevated relative to PAHs in dwarfs. lEngelbracht et alj 
d2005l) also examined low-metallicity systems and found that 
8.0pm emission decreases abruptly relative to 24pm dust emis- 
sion when the metallicity is less than one-third to one-fifth solar. 
M33 is a low-luminosity galaxy hosting HII regions with metal- 
licitie s between solar and one-fifth solar and a shallow rad ial gra- 
dient dMagrini et al.ll2007bt iRosolowskv & Simonll2008l) . PAHs 
are not globally underabundant in M33 although there is a faster 
decline of 8//m emission r elative to longer wavelengths beyond 
3.5 kpc dVerlev et alj|2008l) . It is of interest then to examine pos- 
sible metallicity and radial dependences of PAH features in in- 
dividual HII regions. 

In Section 2 we define young stellar cluster samples in M33 
and analyze their IR properties. In Section 3 we study their UV 
luminosities and colors and define the cluster bolometric lumi- 
nosities. The concept of cluster birthline is given in Section 4 
and tested using M33 young clusters. In Section 4 we discuss 
the implications of our results relative to the dust abundance, to 
the bolometric to Ha luminosity ratio, and to the population of 
massive stars in star-forming sites of M33. In the Appendix we 
analyze possible metallicity dependencies using a small sample 
of young clusters for which gas metal abundances are known, 
as well as dependencies of cluster properties on the mass of the 
associated giant molecular cloud. 



2. Multiwavelength observations of star-forming 
sites 

The closest star-forming galaxies offer a unique opportunity to 
resolve individual star-forming regions. In M33, the star forma- 
tion process can be investigated not only by averaging over se- 
lected areas of the disk or by using very bright complexes, but 
also by observing individual star clusters of moderate size. The 
sensitivity of telescopes such as Spitzer in the IR and GALEX 
in the UV allows us to sample clusters over 4 orders of magni- 
tude in brightness: from the second brightest HII region of the 
Local Group (NGC 604) down to clusters of about 1000 L Q . We 
can look in detail at the site properties where clusters of dif- 
ferent masses are born, and at the emission properties of young 
clusters. We will use photometric information on sources across 
the M33 d isk as given by the S pitzer 8pm and 24pm images de- 
scribed in IVerlev etai] {20071) . b y the GALEX FUV a nd NUV 
images from the GALEX Atlas dbil de Paz et al.ll2007l) . and by 
the Ha image of lHoopes & Walterbosl d2000t) . 



2. 1. The main sample 

We define a sample of young stellar clusters and isolated mas- 
sive stars using the IR c atalog of 51 5 sourc es extracted from the 
24pm Spitzer map by Verle v et ail d2007l) . This catalog has a 
completeness limit of 3 x 10 36 erg s _1 (this limiting luminos- 
ity refers to the 90% level of completeness computed by ex- 
trapolation of the observed cumulative luminosity function at 
L24 > 10 37 erg s _1 ). We remind the reader that the photometry 
for these sources has been done with a varying aperture to match 
the source size. From the catalog we select only sources that 
have an Ha counterpart. We search for Ha counterparts to 24pm 
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sources as follows. We first define a sample of Ha sources by 
extracting 413 sources from the Ho- image using the SExtractor 
software wit h no specified sky coord inate list and with vary- 
ing aperture (Berti n~& Arnoutslll996l specific inputs given in 
Verley et al. 2007). For any source in the IR catalog we then 
define a searching radius, R s , equal to the radius of the source 
in the 24//m image (mino r axis if the s ource PSF is elliptical). 
This radius is derived by IVerlev etaf] d2007l) using SExtractor 
in the framework of the "hybrid" method and varies between 12 
and 200 pc. If at least one of the sources in the Ha sample is 
at a distance less than or equal to the searching radius of an IR 
source, then this infrared source is selected and its Ha flux is 
set equal to the sum of the fluxes of all Ha sources within R s . If 
no Ha sources are found within R s , then we attempt to measure 
the Ha flux using the "hybrid" photometric method devised by 
IVerlev et alj (|2007). With the "hybrid" photometric method, we 
search for Ha emission in the circular area of radius R s around 
the positions of IR sources (specifying the source coordinates to 
the SExtractor software). If again an Ha feature is found within 
the searching area, then this source is selected and its Ha flux set 
equal to the "hybrid" photometric flux. If Ha emission is found 
at a distance larger than R s , then the IR source is not selected for 
our sample. The use of the hybrid method allows us to recover 
many faint Ha counterparts because their approximate positions 
are entered into SExtractor. From both methods we have a final 
sample of 355 sources that have been selected from the 24/im 
Spitzer map and have an Ha counterpart. Of these, only 3 have 
a double counterpart in Ha, less than 1%. We estimate the com- 
pleteness limit of the survey to be 3 x 10 35 erg s in Ha (90% 
level of completeness). 

Using the 8/im Spitzer map, we also search for 8//m counter- 
parts to all 24/mi sources. We employ the same method used for 
finding Ha counterparts. We find that 491 of the 515 sources in 
the 24/im catalog have an 8/im counterpart. For the subsample 
of interest in this paper (the 355 IR sources with an Ha coun- 
terpart), 349, or almost all, have a clear detection at 8/mi. To the 
limit of our completeness for 8/zm, estimated to be 10 37 erg s" 1 
(90% level of completeness), we find that 8% of the sources have 
multiple counterparts at 8/im. If the associated PAHs are located 
along the photodissociation regions and HII shells, then more 
than one 8/im bright rim might be detected per source even if 
one stellar cluster is powering the 24/im emission at the center. 

In Figure 1 we show that the source sizes are always much 
smaller than their average separations, so the sources are dis- 
tinct. The scatter around the average source size is small, com- 
parable to the size of the symbols used in the plot. Source sizes 
decrease radially outward and their average separations increase. 
The 8/mi and Ha images are used at the original resolution, 
higher than that of the 24//m image. The use of images smoothed 
to the spatial resolution of the 24/^m image, 6 arcsecs, would 
have more firmly constrained the total flux at 8/im and Ha for 
some cases, but it would also have diluted the fainter sources 
in the midst of high diffuse emission, making them harder to 
detect. The presence of elliptical source PSFs also makes it dif- 
ficult to match the source fluxes at all wavelengths. Also, for the 
main sample, the 8/im and Ha fluxes might be slightly higher 
in crowded fields than what we derive (i.e., we might be miss- 
ing some multiple faint components with luminosities on the or- 
der of our completeness limit). However, as we shall see later, 
this will not affect the main results derived using the isolated 
source sample. This sample, defined more precisely later in this 
Section, consists of round, compact sources without multiple 
counterparts; they are found by visually inspecting the images 
at all wavelengths. 
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Fig. 1. The average distance between sources, d, and the aver- 
age source size, s, are shown as a function of galctocentric ra- 
dius. The dispersions around the average source sizes are com- 
parable to the symbol sizes. Different symbols refer to mea- 
surements at different wavelengths: 24/mi (open triangles), 8/im 
(open squares), Ha (crosses). 



In Figure 2 we show source luminosities, the 24jum-to-Ha ra- 
tio and the 24/im-to-8/im ratio as functions of the galactocentric 
radius R. The 24//m and 8//m luminosities are defined as vL v , 
where L Y is the luminosity per unit frequency at 24 and 8/im, 
i.e. L v = 4nD 2 F v , F v being the IR flux measured and D the 
distance to M33. L24, La, an d Ln a (uncorrected for extinction) 
are in units of erg s . The ratio Z^/^Hcr has a high dispersion 
and is consistent with having no radial dependence. The slope 
of the linear fit in the log(L24 / 'Lho) - R plane is -0.06 + 0.02 
with a Pearson linear correlation coefficient of -0.14. Brighter 
IR sources lie preferentially at smaller radii, but the radial de- 
pendence is not strong enough to be dominant when the ratio 
of 24/im to Ha luminosity is considered. Ha emission appears 
to have a very marginal radial dependence and a higher disper- 
sion than the IR luminosities. If extinction corrections, generally 
higher at smaller radii, were applied to the Ha source luminosi- 
ties, then the sign of the poor correlation found between L,2aILhc 
and R might be reversed. However, we show in the rest of this 
paper that extinction corre ctions are small; M 33 is a galaxy with 
low dust content (see also IVerlev et alj|2008l) . 

The ratio L24 / La has a radial dependence that is weak but un- 
doubtedly exists: the slope of the linear fit shown in Figure 2 is 
0.10 + 0.01 and the Pearson linear correlation coefficient is 0.34. 
The radial dependence of the 8/im source luminosity is stronger 
than that of the 24/im luminosity, so the 24-to-8/im luminosity 
ratio increases at large galactocentric radii. Because PAHs are 
likely responsible for the 8/mi emission in M33, it might be that 
in the outer disk the lower ISM pressure pushes the photodis- 
sociation regions, where PAHs reside, further from the bulk of 
the 24/im emission. Then our algorithm for finding 8/im coun- 
terparts to 24/im sources wou ld fail at large galactocentric radii. 
However. IVerlev et alj (|2008) have found that the average radial 
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Fig. 2. The ratio Li<\jLn a , La^jL^ and the Ha, 24/vm and 8/vm 
luminosities (in units of erg s ) for IR selected sources in M33 
are shown as functions of their distances from the galaxy center. 
The dashed lines in the upper 3 panels show the survey com- 
pleteness limits. The continuous lines in the bottom 2 panels are 
the best linear fits to the data. 

profile at 8yum in M33 falls off more steeply than that at 24, 70 
and 160/mi, and at UV wavelengths. Hence PAH carriers seem 
effectively underabundant at large radii. There have been claims 
in other galaxies that this might be due to the low metal content 
of the outermost regions, but it is unclear if this a pplies to M33 
becaus e of its very shallow metallicity gradient (Magrini et al. 
l2007bh . 

2.2. Total IR luminosities 

IVerley et al.1 (|2007|) have already pointed out the large scatter be- 
tween the SFR in selected HII regions of M33 computed via the 
Ha optical recombination line and that inferred from the 24/im 
infrared flux. We shall show that this cannot be from extinction 
affecting the Ha line because M33 has a globally low dust con- 
tent. Nor does it result from total IR luminosities of star-forming 
sites i n M33 tha t are no t well represented by 24//m luminosities. 

IVerlev et al.l d2008l) have shown that the total infrared lu- 
minosity (hereafter TIR, defined as the emission between 3- 
1 100 /urn) in star-forming regions correlates with the 24yum emis- 
sion, with a weak additional dependence on the 8-to-24//m ratio. 
They show that for a sample of sources selected at 160^/m, the 
log (F u/Ftir) varies between -0.8 and -1.2 when log (fVi^) 
varies between -0.5 and 0.3. The best fitting relation is 

logL™ =logL 24 + 1.08 + 0.51 logf^^J (1) 

where the 24yum luminosity is vL v , as usual. We shall therefore 
use the above equation to compute the TIR luminosities of star- 
forming sites in our main sample. These are shown in Fig. [3] For 
the 6 sources with no clear 8yum detection, we neglect the last 
term in Eq.Q] In Figure [4] we show the ratio of the TIR to Ha 




Fig. 3. The total infrared luminosity of sources in the main sam- 
ple as a function of the 24/mi luminosity derived according to 
EqH Circles mark sources which belong to the isolated sample. 



flux as a function of the TIR luminosity for all 355 sources in 
the main sample. We encircle sources in the isolated sample (de- 
fined in the next subsection). The sources span over two orders of 
magnitude in the TIR to Ha luminosity ratio and about 2.5 orders 
of magnitude in Ltir. In Figure|4]the Ha flux has been corrected 
for extinction using an average Ay of 1 mag (0.83 mag at Ha). 
Since the maximum extinction we measure in the next Section 
corresponds to Ay =1.7 mag, which is also the average extinc- 
tion m easured towards bright HII regions (llsrael & Kennicuttl 
1980), more heavily extincted than faint ones, we consider a dis- 
persion of 0.25 mag in Ay in the error propagation analysis for 
Figure|4] Extinction uncertainties dominate over photometric er- 
rors for bright sources. 

The TIR-to-Ha luminosity ratio for our main sample has 
a large dispersion and shows no dependence on galactocentric 
radius or source brightness. In the next Sextion we will use 
UV photometry to compute Ha extinction corrections for each 
source in the round and isolated samples. The results prove that 
the large scatter we observe in Figure|4]is intrinsic to the sources 
and is not the result of measurement errors. 

We shall show in the rest of the paper that the scatter in 
Figure [4] is related to three effects: (a) a varying dust opacity 
around sources, which implies that the TIR is not proportional 
to the cluster bolometric luminosity; (b) a non linear relation be- 
tween the bolometric and Ha luminosity of newly born clusters 
when the cluster luminosity is below a certain threshold, and (c) 
cluster aging. 

2.3. UV photometry: the round sample and the isolated 
sample 

We now complement the infrared and Ha photometry with the 
far- and near-UV data. UV colors can also be used to infer the 
age of 24 ^m selected sources. G ALEX-UV images of M33 are 
available dGil de Paz et al.ll2007l) and we shall use them at the 
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Fig. 4. Ratio of the TIR to Ha luminosity as a function of 
the TIR luminosity for IR selected sources. An average value 
Ay — 1 ± 0.25 mag has been used to correct the Ha luminos- 
ity for extinction. The TIR luminosity has been inferred using 
Eq03 Circles mark sources which belong to the isolated sample. 
Dashed lines mark the completeness limits of our survey. 



original resolution (1.5 arcsec) to measure the FUV and NUV lu- 
minosities of 24 pm sources with Ha counterparts. Because we 
perform the UV photometry using a circular aperture, we select 
106 sources at 24 //m using the requirement that the PSF elliptic - 
ity is less than 0.2. We shall call this the round sample. We set the 
radius of the UV photometric aperture equal to the source size 
at 24 pm and consider an annulus with a width equal to 2 pixels 
(3 arcsec for GALEX) around sources for background removal. 
It is not desirable to perform individual UV source extraction 
because of crowding due to the high densit y of UV sources. A s 
pointed out already in previous papers (e.g. Cal zetti et al.l l2005). 
the UV emission peaks are sometimes displaced with respect 
to the infrared emission peaks. Stellar complexes can be made 
of 2 or more populations of different ages, both visible in the 
UV, with the older ones more extended than the young ones. 
In this case the old ones should be included in the background. 
We would like to measure only the UV flux coming from the 
same area covered by IR emission. To strengthen our results we 
have selected a clean sample, made of 26 sources that have a 
clean annulus in both the FUV and NUV images, i.e. no source 
contamination for background subtraction. We shall call this the 
isolated sample. Briefly, the sources in the isolated sample are 
sources from the main sample which satisfy the following crite- 
ria: (0 the PSF of the source at 24pm is nearly circular i.e. the 
ellipticity is < 0.2 and all of the UV flux is inside R s ; (if) the 
source has one Ha counterpart (as defined earlier in this paper); 
(Hi) in the far- and near-UV maps, an annulus of 2 pixels around 
the emitting source is devoid of any sources. 

We will use the round sample, which includes the isolated 
sample, in the rest of the paper to understand some basic prop- 
erties of young stellar clusters in M33. To check that we are not 
oversubtracting the background in the case of sources which do 



not have a clean annulus, we compute the average UV luminosi- 
ties for the isolated sample and for the ensemble of sources in 
the round sample which are not in the isolated sample. The ra- 
tio between the two average luminosities is close to unity, being 
38.556 and 38.365 the FUV and NUV logarithmically averaged 
luminosities (in erg s ) in the round and non-isolated sample, 
and 38.563 and 38.377 respectively for the isolated sample. We 
are therefore confident in the ability of our technique for UV 
background removal to recover the UV luminosities associated 
with young clusters. 

The UV AB-magnitudes were converted to luminosities us- 



ing: 



^uv 



vF v AnD 



(2) 



with D the distance to M33, v = 1.9510 15 Hz for the FUV band, 
v = 1.310 15 Hz for the NUV band, and 



Jy 



10 



23-(AB+48.6)/2.5 



(3) 



Considering that the GALEX flux calibration has been 
revised after the early data release, we apply the recom- 
mended 0.12 magnitude o ffset to the FUV-NUV color defined as 
log [L A (FUV)/L A (NUV)] dBianchi & Efremovall2006l) . We con- 
sider also the small extinction correction, E(B- V) mw - 0.04, 
from the Milky Way in the direction of M33 (Schlege l et alj 
1998). This implies A uv 0.3 for all sources in M33. 



3. UV colors, opacity and bolometric luminosities of 
star-forming regions 

The TIR luminosity of Hll regions is radiation from young stars 
absorbed and re-emitted by dust grains. The leftover cluster ra- 
diation, not absorbed by grains, is mostly emitted in the UV. We 
now infer the extinction of the radiation coming from the gas 
and the stars in the HII reg i ons us ing the observed infrared-to- 
far UV ratio. As in lCalzettD d200ll) . we shall use the factor 1.68 
to account for the bolometric correction of the stellar emission 
relative to the FUV. 

Since our IR and UV data refer to clusters still embedded in 
HII regions, the dust optical depth for the stellar continuum radi- 
ation is similar to that for the Ha line emission. We are not con- 
sidering the UV light from previous generations of stars. We in- 
fer the FUV extinctions after measuring the TIR and FUV source 
luminosities, and then we relate this to the extinction at other 
wavelengths as in ICalzettil d200lh : 



5.7 x log 



1 L 



TIR 



1.68 L r 



+ 1 



i-NUV 



= 0.72xA 



FUV 



Ay = 0.3 xA, 



A Ha = 0.25 xA f 



(4) 



(5) 
(6) 
(7) 



We show in Figure|5]the ratio LjirILfuv versus the GALEX 
color L / i(FUV)/L l i(NUV) for sources in the rou nd sample. By 
comp aring this Figure with Fig. 9 and 10 of ICalzetti et al.l 
(120051) . we see that our selected sources have all log (L T[R /L FUV ) 
between -1 and 1.2, corresponding to Ay < 1.7. The derived 
value of Ay increases on average towards bright sources. The 
average UV color is on the order 0.3 mag. All sources in the 
isolated sample and nearly all sources in the round sample have 
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Fig. 5. Opacity expressed as infrared-to-FUV luminosity ratio 
versus the observed UV color for all sources in the round sample. 
Circles indicates sources in the isolated sample. All sources to 
the left of the line marked 10 Myrs are younger than 10 Myrs 
according to the instantaneous burst model. All sources below 
the mark Ay = 1 .7 have visual extinction Ay < 1 .7 mag. 

UV colors > 0.2 mag., i.e. they are younger than 10 Myrs ac- 
cording to the instantaneous burst model for sola r metallicity 
and for a starburst dust distribution (Calzetti et al. 2000). Exact 
ages will depend on the metallicity, burst model and dust dis- 
tribution. An extinction curve like the one observed in the Small 
Magellanic Cloud, as well as the simultaneous presence of newly 
formed massive stars and an older UV-emitting population, im- 
ply even younger ages for Ha emitting sources in our selected 
star-forming regions. 

We estimate the bolometric luminosity of our sources as: 

Ltd = L T1R + [v FUV L v (FUV) + v NUV L v (NUV)] (8) 

where vfuv and v^uv are the effective frequencies of the 
GALEX bands and the UV luminosities are uncorrected for ex- 
tinction. The bracketed quantity estimates the cluster UV lu- 
minosity, between 912 and 3000 A. In the absence of dust the 
UV luminosity is a good approximation of the bolometric lumi- 
nosity for clusters with ages betw een 1-10 Myr, indep endent of 
their mass (as from Star burst99 Leith erer et al.l ll999) (see also 
iBell et aLll2005t iThilker e t al. 2007] for similar methods used to 
estimate Lboi)- We shall give another proof in the next Section 
that indeed our selected clusters fall into this age category. When 
dust absorb part of the UV radiation and re-irradiate it in the IR, 
Ltir is needed in the right hand side of the above equation to 
correctly estimate Lj /. 

We apply extinction corrections to the FUV and Ha lumi- 
nosities of sources in the round sample according to the model 
described earlier in this Section. The good correlation between 
extinction corrected UV luminosities and 8//m luminosities is 
evident in Figure [6] This is expected since PAHs are excited by 
the UV radiation of the nearest star-forming region. This nice 
correlation confirms the goodness of our extinction corrections. 



Fig. 6. Ha and 8/im luminosities are shown as functions of the 
FUV luminosities for all sources in the round sample. Circles 
indicate sources in the isolated sample. Luminosities are in units 
of erg s _1 ; Ln a and Ls have been corrected for extinction using 
the source TIR-to-FUV luminosity ratio. 

Thil ker et alj {2007) have found that PAH emission at 8//m is 
suppressed within strong star-forming regions. The linear cor- 
relation between extinction corrected far-UV luminosities and 
8/mi luminosities, which does not break at high luminosities, 
suggests that this is not happening in M33. 

A correlation between FUV and Ha luminosities is also 
present. Given the different dependencies of these two luminosi- 
ties on cluster mass and age, the presence of some scatter in their 
relation is not surprising. The next section will better address this 
issue. 

4. The cluster birthline: from bright to dim candles 

Because we resolve stars and clusters in M33 at IR and optical 
wavelengths when the brightest member is earlier than a B2 type 
star, the expected relation between IR and Ha emission is far 
from linear even under the circumstance that most of the bolo- 
metric luminosity of massive stars is absorbed by dust and re- 
emitted in the IR. This is a result of the different depend encies of 
the bolometric and Ha luminosities on the stel lar mass {Panagia 
119731: IVacca et alj|1996t ISternberg et alj|2003l) . considering the 
lack of the highest mass stars in some clusters. In this Section 
we would like to investigate in detail the emission by individual 
clusters at different wavelengths to judge the reliability of using 
infrared and Ha emission as a quantitative tracer of recent star 
formation for a range of luminosities. 

4. 1 . Modelling cluster birthlines 

To compute the Ha luminosity of a stellar cluster, one needs to 
know the mass distribution at the upper end of the IMF. Even 
though the stellar IMF is well known for inte rmediate stellar 
masses, where the original exponent derived by Salp eterl d 1 95 5b 
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is now widely tested and used, uncertainties remain at the highl- 
and low-mass ends. 

The IMF may be interpreted as a probability distri- 
bution function that specifies the probability for a ran- 
domly chosen star to lie within a certain mass range (e.g. 
ICervino & Lur idiana 2006; Barker et al.ll2.Q08l for stochasticity 
in the IMF). The absence of stars with very high masses sug- 
gests there is a fundamental upper limit to the IMF (> 120 
150 M see Massey & Huntei 



and this star will be of mass M m . If £(M) is the number of stars 
produced per unit mass, 



ip per limit to the IMt (> 1ZU - 
r 19983: lElmegreenl 120001 iFigei 



2005: lOev & Clarke - 2005: Koen 2006; Hoversten & Glazebrook 
2008). The most massive clusters sample the IMF out to this 
maximum possible stellar mass and no further, even though there 
are enough stars to do so if the IMF continued with the same 
slope. Low and intermediate mass clusters do not generally have 
stars near the maximum possible mass, but their most massive 
stars tend to h ave a mass that increases with the cluster mass 
dLarsonl ll982). This trend leads to the question of whether the 
upper limit to the stellar mass in a particular cluster depends 
on the cluster mass because of some physically limiting pro- 
cess, or whether it follows only from randomly sampling the 
IMF. In the first case, low mass clusters could not make high 
mass stars. In the second case they could, as long as there is 
enough gas, and intermediate-mass clusters should occasionally 
be found with unusually massive stars - "outliers" in the IMF. 
An important difference arises for the summed IMF from many 
clusters: it should be steeper than any cluster IMF in the first 
case and the same as the average cluster IMF in the second case. 
The summed IMF is the basis for the present day mass function 
in whole galaxies and is therefore important to understand for 
population models and studies of galaxy evolution. A recent de- 
bate in the literature suggests the question a bout maximum stel- 
lar mass should be settled ob servationally (We idner & R roupa 
120041 120061: lElmegreedl2006l) . 

In order to model the distribution of stellar masses within 
a young cluster, we follow two different approaches using the 
same IMF. In one case we assume that each cluster has a max- 
imum stellar mass that is explicitly related to the cluster mass. 
We call this the maximum mass case. In the second case, we as- 
sume that stars of all masses can form in clusters of all masses, 
provided there is enough gas, and that the choice of a particular 
stellar mass is random, following the IMF. We call this the ran- 
domly sampled case. Once we specify how a cluster is built we 
need to establish L c } , and , the cluster bolometric and Ha lu- 
minosity. This is done by summing the contributions of all stars 
in the cluster using stellar masses, Ly-continu um, and bolometric 
luminosities tabulated b v lVacca et alj (1 19961) for stellar masses 
87.6<M<19.3 M Q and by|Panagia|(ll973|) for 19.3<M<8 M G . We 
extrapolate luminosity values logarithmically if stellar masses 
are outside these mass ranges and interpolate logarithmically 
for masses that are intermediate to those shown in their tables. 
We assume ZAMS (Zero Age Main Sequence) and solar metal- 
licity and divide the number of Ly-continuum photon rates by 
7.3 x 10 11 ph erg -1 to obtain Ha luminosities. 

In what follows we compute the theoretical cluster birthline 
in the plane \og(Lboil LH a )-\og{Lb i) using a stellar IMF with a 
Salpeter slope a = -2.35 between 1 and 120 M Q . No stars with 
mass lower than 1 M Q or higher than 120 M will be considered. 
The IMF in stellar clusters extends down to masses lower than 
1 M Q but these stars give a negligible contribution to the cluster 
luminosity, which is the relevant quantity for this Paper. 

In the maximum mass case, given a number N of cluster 
members with masses above 1 M Q , there will be a certain mass 
limit Mm above which the probability of finding one star is unity 



%{M)dM = CM a dM, 



(9) 



then the number of cluster members with masses above 1 M , 
and the cluster mass and luminosity contributed by those stars, 
can be written as 



» C MN M a dM 
N;/ = I CM°dM = 1 + J1 



rl20 . , 
JMu 



CM a+l dM = M, + 



^ N M a+l dM 

rl20 ~ ,„ , 
JM N 



r \20 



CL(M)M a dL = L(M t ) 



JJ * L(M)M a dL 
[! A 2 °M a dM 



(10) 



(11) 



(12) 



We computed the normalization constant C by assuming that the 
number of stars between 120 M and Mjv is unity and we con- 
sidered the IMF from M N to 1 M as a continuous d istribution 
function with or = -2.35 (e.g. lWeidner & Kroupal2006l and ref- 
erences therein). We computed the cluster masses and luminosi- 
ties in this way for a sequence of clusters with maximum stel- 
lar masses corresponding to ZAMS stars of spectral types B2, 
Bl, .... 04, and 03. We set Mjv equal to the mass of these stars 
and consider also that M», the mass of the star born in the mass 
range Mn-120 M , is equal to the lower limit of this range, M#. 
Statistically M^ is the most likely value for mass in the interval 
Mat-120 M , but the assumption that for a given cluster mass 
no stars are ever made with M>M^ is stronger. In Figure [7] we 
plot with filled circles the cluster birthline made with the max- 
imum mass assumption. That is, each filled circle corresponds 
to a spectral type in the above sequence, with the corresponding 
bolometric-to-Ha- cluster luminosity ratio on the ordinate and 
the bolometric luminosity on the abscissa for clusters that have 
a maximum stellar mass with that spectral type. 

In the randomly sampled case, there is no explicit relation- 
ship between cluster mass and maximum stellar mass as in equa- 
tions (TlOt - CfT2b . Even though M^ is the most likely mass for the 
cluster's most-massive member, the actual mass can be as large 
as 120 M , although the probability that such a massive star ac- 
tually forms decreases as the cluster mass decreases. To model 
this case, we simulate 40000 clusters that are distributed in num- 
ber according to their mass between 20 and 10000 M as 



N( M cl ) dM d = D cl M'l, dM c/ 



(13) 



where D c i is a constant and 6 is th e spectral index of the ICMF, 
usually assumed to be 6 — —2 (e.g. lde Griis & Parmen tier 2007, 
and references therein). We then populate each cluster with stars, 
randomly selected from the stellar IMF described above, un- 
til the simulated cluster mass gets above the selected cluster 
mass. We compute the bolometric and Ha luminosity of each 
cluster using the same ZAMS models as in the maximum mass 
case. Figure [8] shows the resulting distribution of L^/L^o. ver- 
sus Lbd with one plotted point for each modelled cluster. The 
lower boundary of the distribution is where single massive stars 
lie. Our simulation takes into account incomplete clusters, such 
as single OB associations born in the field or single stars scat- 
tered out of loose clusters. Notice the presence of a second edge 
of high density points in the distribution. This is from a change 
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Fig. 7. Filled circles mark the expected values of the bolometric - 
to-Ho' luminosity ratio for the maximum mass case. Open and 
filled squares are the mean and median values respectively of the 
same ratio for the randomly sampled case using a cluster mass 
function index of 6 — -2. Errorbars point out the dispersions 
around the mean and median values. 



in slope of the stellar L# a -L£ / relation, which produces a cusp 
in the diagram. In bins of equal size in log L/,,,/ units, we plot in 
Figure |7] the average and the median of the distribution for the 
randomly sampled case and we plot their 1-cr values (enclosing 
16.5% and 83.5% of the points in each bin). Notice that the ran- 
domly sampled case predicts median values and average values 
of log Lbo/fLna that are below those of the maximum mass case. 
This is a result of the occasional presence of bright stars or out- 
liers even in low-mass clusters. 

If low mass stars form first, a cluster enters the birthline from 
the top left and follows it down and to the right as L/, a / increases. 
If, instead, massive stars are born at random times during clus- 
ter formation, then the cluster can jump down on the simulated 
birthline and then move towards slightly higher values of L^,/ 
and LboifLfja- This happens until all of the gas available to fuel 
star formation is used or SF is quenched by feedback. Cluster 
aging from the death of massive stars increases the Lfo,//L// ff ra- 
tio and then the cluster moves upward from the birthline. This 
is because the death of massive stars makes the cluster Ha lu- 
minosity fade away more rapidly than bolometric luminosity. 
Using the public code Starburst 99, we evaluate how the clus- 
ters evolve off the birthline. The vertical dashed lines in Figure 
[8] show the aging effect. A cluster reaches the top of the Figure 
in about 10 Myrs. Two fundamental properties of the predicted 
cluster birthline are: (a) the linear correlation between L& / and 
Lua, which holds for very young massive clusters (the horizontal 
part of the birthline), breaks down below a critical value of the 
cluster bolometric luminosity, depending on the stellar IMF at 
the high mass-end; (b) all stellar clusters should lie on or above 
the birthline shown in Figure|7]for the two cases examined. 




37 38 39 40 41 

Fig. 8. The distribution of the ratio of the bolometric-to-Ho' lu- 
minosity versus the bolometric luminosity for newly born clus- 
ters simulated in the randomly sampled case. Vertical dashed 
lines mark evolutionary sequences, i.e., aging moves clusters up- 
ward from the simulated birthline because Ha luminosities fade 
away more rapidly than bolometric luminosities. Clusters reach 
the top of the plot (Lhoil^Ha — 4) in about 10 Myrs. The dash- 
dotted lines are indicative of when the most massive stars in a 
cluster leave the main sequence. 



4.2. The M33 test case for the birthline 

If most of the luminosity of young clusters were absorbed by 
dust and re-irradiated in the IR, then the Ltir - L^ / data points 
in the log(L r/ ^/L Hu ,) - log(L r/R ) plane should all lie above the 
birthline sequence. However, by looking at Figure|9] which plots 
the observed values on the birthline plane, we see a clear discrep- 
ancy with the birthline theoretical prediction. Many clusters lie 
below the birthline, for both the maximum mass case and the 
randomly sampled case. It has to be pointed out that any heav- 
ier extinction correction, as well as corrections due to the loss 
of ionizing photons leaking out of the HII regions, will move 
the data points further down, toward lower values of L;, //L//q,. 
Since any data point below the birthline is unexplained we must 
conclude that our assumption, namely that the TIR luminosity 
traces the bolometric luminosity, is not correct for most of the 
clusters observed. This discrepancy with the birthline cannot be 
the result of an unreliable 24pm - TIR conversion because the 
dispersion in the 24-pm to TIR luminosity ratio is only about 0. 1 
in the log. Also recall that individual extinction corrections have 
been applied to Ha luminosities (see Section 3). 

As outlined in the previous Section, we now complement IR 
with UV photometry to derive the cluster bolometric luminosity. 
In Figure[10]we plot the L\, i-\a-Lji a ratio as a function of Li,„i for 
the round sample marking (with circles) members of the isolated 
sample. Clusters in the round and isolated samples lie above the 
birthline predicted by the randomly sampled case with a very 
few exceptions. Members of the isolated sample, for which UV 
photometry is more accurate, follow closely the birthline and 
this underlines their young age. Only one cluster in the isolated 
sample with L& / - 10 40 erg s _1 seems to lie below the birthline, 
having a lower value of L£ //L# ff . The most likely explanation 
of the Ha excess is that this cluster has formed a star more mas- 
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Fig. 9. The total infrared-to-Ha luminosity ratio for data in the 
round sample and in the cleanest isolated sample (symbols sur- 
rounded by circles). Under the assumption that Ltir=Li, i, we 
can compare the cluster data with the predicted cluster birth- 
lines for the randomly sampled case (filled squares) and for the 
maximum mass case (large filled circles). Individual extinction 
corrections have been applied to Ha luminosities (see Section 
3). The selected samples of young clusters is clearly not along 
or above any birthline, as it should be. Hence the assumption 
LriR-Ltoi is incorrect. 




37 38 39 40 41 

log L bol [erg s- 1 ] 

Fig. 10. Bolometric to Her luminosity ratio versus bolometric lu- 
mosity for data in the round sample and in the cleanest isolated 
sample (dots surrounded by circles). Bolometric luminosities 
have been estimated according to Eq.|8]and individual extinction 
corrections have been applied to Ha luminosities (see Section 3). 
The predicted cluster birthlines for the randomly sampled case 
(filled squares) and for the maximum mass case (filled circles) 
are also shown. The dash-dotted line indicates the cluster posi- 
tion when the most massive stars for the randomly sampled case 
leave the main sequence. 



sive than 120 M Q , but additional analysis is needed before draw- 
ing any definitive conclusion on this. In general the agreement 
is good also for the round sample: almost all of the clusters lie 
above the birthline for the randomly sampled case. This implies 
that the round sample is made up of relatively young clusters 
too. In Figure [10] we add the main sequence boundary (dash- 
dotted line). This marks the value of L& //L#q, f° r a cluster at the 
end of the main sequ ence lifetime for single massive stars (e.g. 
IStahler & Pallal2005l) . It is impressive how the isolated sample is 
bounded by the birthline at the bottom and by the main sequence 
boundary at the top. That means that effectively we have traced 
the sequence of young compact clusters before supernovae dis- 
rupt them and blow away dust grains. 

If the birth of stars in clusters is regulated by the statistical 
character of the IMF and there is no physical or absolute relation 
between the maximum stellar mass in a cluster and the cluster 
mass, then the predicted cluster birthline is in good agreement 
with the data. This agreement is illustrated by the relative num- 
ber of sources below the birthline, where there should be none. 
Only 5% of the clusters are more than 3-cr below the birthline 
in the randomly-sampled case, whereas 25% of the clusters are 
more than 3-cr below the birthline in the maximum-mass case. 
This result favours the randomly sampled model for the stellar 
population in clusters. 

We infer that some of the clusters with the lowest Li, //L# (! , 
ratios for their L/, / values have outlier massive stars, i.e. massive 
stars without their usual proportion of low mass stars. In these 
cases, a star has formed with a mass that is larger than the av- 



erage upper mass limit for a cluster of that mass. It would be 
interesting to study these cases more to see if there are other un- 
usual characteristics, such as extreme mass segregation, cloud 
disruption, or cloud temperatures. 

The simulated hboil^Ha distribution is very marginally de- 
pendent on the value of S for the cluster mass function. For 
-1 < 5 < -3 the distribution of the median value is within 
the 1-cr error bar value s of the more widely used 5 = -2 distri- 
bution (see FigurefTTI i dHunter et alj|2003l Ichandar et al.| [2006: 
iDowell et alj|2008l) . Figure [12] shows the effect of using an IMF 
index steeper or flatter than the classical Salpeter value. There 
is a better agreement of some cluster data with a steeper IMF at 
the high-mass end, but the number of clusters below the birthline 
increases. So, we favor a Salpeter IMF with some cluster aging 
to bring them above the birthline. A diagram of L^/ versus UV 
colors shows that the most luminous sources are also the oldest 
in our isolated sample. This accounts for the slight rise in plotted 
points above the birthline for large L^/ in Figure[10] Additional 
data is needed to better constrain the ages of those clusters. 

5. Summary and Discussion 

M33 is ideal for deriving the properties of young star-forming 
clusters and their massive stellar populations. It contains star- 
forming complexes with a wide range of luminosities, and it is 
close enough that even clusters with one O- or B-type star can be 
localized. We have examined in this paper young stellar clusters 
selected in the 24yimi map of M33 that have Ha counterparts. The 
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Fig. 11. Median values for the randomly sampled case obtained 
using different indices for the cluster distribution function: open 
squares, filled squares and open triangles correspond to 6 — 
— 1, -2, -3 respectively. Filled dots show the data for the round 
sample. Data for the cleanest isolated sample are surrounded by 
circles. 
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Fig. 12. Median values for the randomly sampled case obtained 
using different index for the stellar IMF: filled triangles, filled 
squares and open square correspond to a = -2.80, -2.35, -1 .50. 
Filled dots show the data for the round sample. Data for the 
cleanest isolated sample are surrounded by circles. 



8-to-24yum luminosity ratio has a small dispersion and shows a 
dependence on the galactocentric radius. The total IR-to-Ha lu- 
minosity ratio, Ltir/Lhu, shows a larger dispersion, especially 
towards faint sources, and no correlation with galactocentric dis- 
tance or source IR luminosity. The large scatter in the TIR-to-Ha 
flux arises mostly because of variations in the local dust abun- 
dance which adds to some aging effect. Only a certain fraction of 
the bolometric luminosity of 24/im-selected sources can be ab- 
sorbed by grains and re-emitted locally at IR wavelengths. This 
fraction might be small for low-luminosity sources that are not 
born along the spiral arms of M33. 

We have calculated the bolometric luminosity for 106 
"round" stellar clusters, L 00 u using UV and IR photometry. This 
has allowed us to test the concept of the cluster birthline intro- 
duced in this paper. The cluster birthline, defined in the param- 
eter space log(L Ao ;)~log(Lfc //Lffo), is the line of birth of young 
clusters, a theoretical lower boundary for the ratio L BO ifLji a for 
each Lboi. We show that along the cluster birthline the relation 
between L& ; and L# a is linear for high-L& ; clusters but is non- 
linear for clusters with L DO i < 3 x 10 39 erg s _1 . Deviations from 
linearity are expected because of two effects: the IMF is not fully 
sampled in low mass clusters, and the number of photons emitted 
by a star with energies above the hydrogen ionization threshold 
has a stronger dependence o n stellar mass than does the stel - 
lar bolometric luminosity (e.g. Panagi afl973l : IVacca et alj l996). 
For high L DO i clusters, the IMF is usually sampled out to the 
highest possible stellar mass and then the ratio of high-to-low 
mass stars is constant, making the L# ff to L DO i ratio constant. 
For low Lboi clusters, the IMF becomes depleted of the high- 
est mass stars; these stars have the lowest ratio of Lfc //L# a and 
their depletion raises the average ratio. The observed flatness of 
LboifLfja for luminous sources implies either that there is a maxi- 



mum possible stellar mass equal to ~ 120 M Q , as assumed in our 
models, or that the L# ff -to-L6 / ratio for individual stars becomes 
constant at M > 120 M Q . 

We modeled the cluster birthline by populating clusters of 
various masses with stars having masses selected from an IMF. 
Several cluster mass function slopes and IMF slopes were used 
for comparison. We modeled the upper end of the IMF in two 
cases. The maximum mass case assumed that each cluster pro- 
duced stars with masses only up to a maximum value that equals 
the average maximum stellar mass for a cluster of that mass. In 
this case, low mass clusters can produce only low mass stars. 
The randomly sampled case assumed that each cluster can pro- 
duce stars over the full range of the IMF, in which case a low 
mass cluster can occasionally produce a high mass star (pro- 
vided the cluster mass exceeds the stellar mass, of course). Both 
cases have the same average maximum stellar mass for clusters 
of each mass, but only the second case can produce outlier stars 
- stars with masses significantly above the average maximum for 
that cluster, which means above the upper end of the smooth and 
declining part of the IMF that is observed at intermediate stellar 
mass. Comparing the two cases, the clusters birthline in the max- 
imum mass case has a higher value of L BO i!Ln a over the range of 
Lboi where O- and B-type stars are starting to populate the upper 
end of the IMF. This range extends from Lboi slightly less than 
10 38 erg s _1 to Lboi slightly larger than 10 40 erg s . At smaller 
Lboi, stars producing H<? emission do not form in the maximum 
mass case and are highly unlikely in the randomly sampled case. 
At larger Lboi, me IMF is fully sampled in both cases. 

From the round sample we selected 26 sources which appear 
isolated and compact in the UV maps i.e. the annular region 
around them, used for background subtraction, is not contami- 
nated by UV sources. All but one of the stellar clusters in this 
sample, called the isolated sample are compatible with cluster 
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birthline predicted in the randomly sampled case for the Salpeter 
IMF with maximum stellar mass limit of 120 M Q . For the iso- 
lated sample in fact 25 out of 26 clusters lie on or above this 
birthline. Considering the whole round sample we find similar 
fractions: 100 over 106 clusters are compatible with the same 
birthline. In contrast, 25% of the clusters in the round sample 
lie below the theoretical birthline in the maximum mass case, 
and hence are incompatible with this cluster population model. 
M33 has therefore provided a positive test to the cluster birth- 
line concept introduced in this paper. Stellar clusters are born 
along the birthline and aging moves the clusters above it. The ob- 
servations also suggest that stars randomly sample the IMF and 
that clusters occasionally have outlier massive stars. Clusters in 
the isolated sample lie below the upper main sequence bound- 
ary (t < 3 - 4 Myr), which is consistent with their compactness 
and young age. Below this boundary, the most massive stars that 
a cluster formed are still present (assuming that all stars form 
at about the same time), so L^z/L/^ remains low, and super- 
novae have not yet disrupted the cluster or blown away the dust. 
No cluster in the isolated sample lies above the upper main se- 
quence boundary, but some clusters in the more numerous round 
sample do. This difference might indicate that some members of 
the round sample have lost their most massive stars. However, in 
the opacity-UV-color diagram, the two samples occupy the same 
area so their ages are about the same. An alternative explanation 
for the high L^/L/^ ratio is then leakage of ionizing photons 
from HII regions. Clusters in the round sample are in fact not as 
compact and isolated as the clusters in the isolated sample, so 
the non-isolated clusters could be more prone to leakage. Future 
studies that resolve individual stars should help to answer this 
age-versus-leakage question. If massive stars are observed di- 
rectly, and in proportional to the IMF, then leakage must be the 
cause of the higher Lf,„//L# ff ratio. 

Luminosities at 8 and 24yum correlate with gas metal abun- 
dances; the dispersion is large because there are few sources with 
an accurate metal abundance determination. While we did not 
find a significative sample of 24/im sources embedded in GMCs 
which are not detected in Ha, we find 34 sources associated with 
GMCs which have an Ha counterpart. The 24yum luminosities 
of these sources correlate with associated GMC masses, even 
though with a large scatter. 

From studies in the solar neighborhood we know that a stel- 
lar cluster in its youngest stage is likely to be deeply embed- 
ded in molecular clouds. Later, mechanical and radiative effects 
of the most massive stars disrupt the cloud, the extinction de- 
creases and the cluster becomes visible in the ultraviolet. The 
time interval in which massive stars are detectable through IR, 
Ha and UV emission depends not only on stellar lifetimes but 
also on the opacity of the parent cloud and on its time evo- 
lution. The characteristics of molecular clouds seem linked to 
the large-scale morphology of the galaxy. Earlier studies of a 
radio-selected sample of thermal sources in M33 have provided 
1 1 you ng, optically-visible ste llar clusters but no embedded star 
cluster (Buck alew et alj|2006l) . A small sample of compact in- 
frared selected sources with no Ha counterpart has been ob- 
served in the CO J=l-0 and J=2-l lines (Corbelli et al. 2009, 
in preparation). The weakness or absence of CO lines suggests 
that these sources are not embedded proto-stellar clusters in the 
process of formation but a more evolved population. In M33 the 
amount of extinction seems generally lower than in the Milky 
Way. The absence of large molecular complexes and a steeper 
mass spectrum implies that even in an early SF phase many stel- 
lar complexes may not be highly obscured. Later, winds and SN 
explosions remove efficiently the dusty envelope and the cluster 



fades away in the IR before less massive stars get off the main 
sequence. 

The concept of a cluster birthline together with the high res- 
olution of future telescopes seems a promising way to analyze 
star formation in external galaxies. 
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Appendix A: The metallicity and GMC samples 

We now investigate if properties of star-forming sites in M33 
are related to the host gas metallicity or to the parent molecu- 
lar cloud mass. We first conside r a sample of HII regions ob- 
served by Ma grini et alj d2007bl) for which metallicities have 
been determined via optical spectroscopy. Extinction can be es- 
timated via the Balmer decrement, even though this is only an 
upper limit to the overall extinction through a whole HII region. 
Unfortunately the sample of HII regions with known metallic- 
ity for which radio continuum or Paschen lines are available for 
determining the extinction deeply into the star-forming region is 
limited to 5 bright sources and therefore the sample is not sta- 
tistically meaningful. Our metallicity sample is made of 3 1 HII 
regions with known O/H abundance, Balmer decrement, and ac- 
curate IR photometry at 8 and 24 /im. Ha emission has been cor- 
rected for extinction using the Balmer decrement and the relation 
Ay — 1.2 x Au a . The extinction (which gives values Ay < 1.7 in 
the whole sample) is comparable to what has been derived with 
more accurate estimates. There is no trend of extinction with 
metallicity or with gas column density that is associated with 
the source site. Low column density gas hosts preferentially low 
metallicity HII regions. 

The sample of sources with IR counterparts shows a shal- 
low radial metallicity gradient (Fig. |A. lb . If we restrict our sam- 
ple to sources where the electron temperature diagnostic lines 
have been detected, and hence the electron temperature can be 
measured (filled dots in Figures |AT| and |A~2l , the n the gr adient 
is compatible with that derived by MagrinietalJ ((2007b) using 
the same selection criteria. The L^aIL^ and Li\IL}i a ratios do 
not show any clear metallicity dependence because sources with 
high metallicity (say O/H above -3.5) have high luminosities at 
optical and IR wavelengths (Fig. |A.2| i. The dependence is how- 
ever very marginal if we consider only sources with electron 
temperature determinations, whose metallicities are more cer- 
tain. The scatter around the mean values increases as the metal- 
licity decreases, but this might be an effect of large galactocen- 
tric radii, where most of the low metallicity regions are found. 

IR luminosities show a dependence on metallicity. Bright 
sources at 8 and 24/im are born where metallicity is high while 
low metallicity gas host a wider range of source luminosities. 
However uncertainties in the metallicities of bright sources that 
have no detections of the electron temperature diagnostic lines 
are large and do not allow us to draw firm conclusions. 

We also analyze the BIMA giant molecul ar cloud dataset 
(13" spatial resolution; Eng argiola et alj 12003) to see if any 
sources in the main sample are associated with GMCs and in 
these cases if the source luminosities are related to the parent 
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Fig.A.l. Metallicity as a function of galactocentric radius for 
sources with IR counterparts. Filled dots represent sources 
where the electron temperature diagnostic lines have been de- 
tected. The metallicity measured in terms of O/H from a com- 
pilation by iMagrini et alj d2007alfbh is also plotted as a func- 
tion of the ratios of the 24/vm to the Ha luminosities, and the 
T Ajiva to the 8;um lumin osities. Error bars in log(0/H) are given 
in Magrin iet al.l (l2007allbl) and are on the order 0. 1 . Error bars in 
the luminosity ratios are smaller than the symbol sizes. 
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Fig. A.2. The metallicity measured in terms of O/H from a com- 
pilation of Magrini et al. (2007a b) is plotted as a function of the 
24jt/m, 8/im, and Ha luminosities for sources with IR counter- 
parts. Filled dots represent sources where the electron tempera- 
ture diagnostic lines have been det ected. Error bars in log(0/H) 
are given in Magri nTet alj (|2007a b) and are on the order 0.1. 
Error bars in the luminosity ratios are smaller than the symbol 
sizes. 



GMC masses. Figure lA~3l shows the GMC masses as functions 
of the L24 and L^a/Lho values for sources within 6.5 arcseconds 
of a GMC. This limiting radius was chosen to m atch the BIMA 
spatial resolution of the lEngargiol a et al.l (2003) GMCs survey 
that was used to identify the clouds. There is a dependence of the 
IR source luminosity on the GMC mass, despite the large scat- 
ter. The slope of the correlation is 0. 17+0.05 with Pearson linear 
correlation coefficient of 0.53. The ratio LqaILhc f° r the GMC 
sample varies by about 1 .5 order of magnitudes for sources asso- 
ciated with small clouds, and by even less for sources associated 
with the most massive clouds. Sources in the GMC sample are 
then confined to a smaller range of L24/Lh than that reported in 
Figure DOl for the full sample. The ratio Z^/Lhq. shows no clear 
dependence on the GMC mass. 
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